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u Aktixct Recent research indicates that the use of remote sensing • 

techniques for the measurement of near surface soil moisture could be I 
practical in the not too distant future. Other research has shov'n that' 
^nfi ItrPiion rates, especially for average or frequent rainfall events, j 
are extremely sensitive to the proper definition and con.sideration 6f j 
the role of the soil moisture at the beginning of the rainfall Thus, j 
dt is important that an easy to use, but theoretically sound, rainfall 1 
iinfiltrat ion model be available if the anticipated remotely sensed | 
soil moisture data is to be optimally utilized for hydrologic Simula- I 
rion. I 

I A series of numerical experiments with the Richards' equation for 
Sn array of conditions anticipated in watershed hydrology were used to I 
develop functional relationships that describe temporal infiltration | 
rates as a function of soil type and initial moisture conditions. ' 
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CHAPTER 1 


INTRODUCTION 

Soil Moisture conditions in the near surface soil regime are importaiit 
drivers to many key hydrologic piocesses such as infiltration, moisture re- 
distribution, and evapo-transpiration. Recent progress in microwave and thermal 
inertia techniques show great promise as tools for defining the spatial as well 
as temporal variations in soil moisture. Host hydrologic models in present use, 
however, were designed without recognizing the growth of remote sensing tech- 
nology and are not structured to take full advantage of this information. In 
fact, models were deliberately simplified in their original development because 
of the absence of spatial and temporal soil moisture information that remote 
sensing technology may be able to provide. 

The potential use of soil moisture measurements in the proper definition 
of watershed rainfall-runoff predictions was demonstrated in a recent study by 
Milkening (1981) using a physically-based mathematical model to study the in- 
filtration process. The results of Nilkening's numerical simulations showed that 
rainfall excess volumes from single event storms are quite sensitive to antecedent 
or initial moisture conditions, especially in more frequent rainfall events. 
Wilkening snowed that rainfall excess predictions were quite sensitive to near 
surface moisture conditions, soil hydraulic properties and storm characteristics 
.'•uch as intensity, duration and temporal distribution. Using such a rigorous 
physically-based model, the impact of initial moisture conditions on the infiltra- 
tion process and resulting rainfall excess can be properly represented. The use 
of this type of model is not practical on a watershed scale, however, because 
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of the exc«>sivc coaputational costs involved In solving the governing physical 
equations an J the large aaount of spatially varying data required to represent 
the response of soils Including the soil aoisture-retentlon curve, the 
saturated conductivity and the unsaturated conductivity function. 

Rather then atteapt to directly Incorporate rigorous, but cuabersome, 
Bathec.allcal aodels into applications oriented hydrologic siaulation, an alternative 
approach is to conduct nuaerical experiaents with such aodels and analyze 
the results to establish relatively siaple functional relationships that 
accurately describe infiltration capacity within soae range of conditions. 

These functional relationships, describing iifiltration capacity in terns of 
initial moisture conditions, general soil texture classifications, and accomodating 
temporal rainfall variations, would be computational ly efficient and could 
significantly improve current infiltration components within both single event 
and continuous simulation models. In particular, these theoretically sound, 
yet siaple infiltration models that logically incorporate' the importance of 
initial soil moisture conditions would be of treoMndons value when interfaced 
with accurate remotrly- sensed soil moisture measurements made on watershed 
basis. 

Several efforts have conducted that use numerical sionilations of theoretical 
infiltration events as a tool to develop practical methods for use in watershed 
analysis. Smith (1972), demonstrated a process of curve fitting to describe the 
properties of a complex theoretical model of porous media flow. The para- 
meters of a parametric equation were fit to time varying infiltration capacity 
curves that resulted from numerical simulations using experimental soil data 
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froB Btvcral difftrmt lolls. Bloo«fl*ld «t al., (1981) used mia«ric«l aiaulations 
to oxaBln* the validity of aapirical equation! relating infiltration capacity 
to toil Boiature itorage in the upper soil zone. The Blooafield group demon- 
strated that serious conceptual errors in traditional Infiltration equations 
could be iaproved when the paraaeters were fit to infiltration-storage relation- 
ships derived froa nuaerical siaulations. 

The objective of the present study was to develop a siaple but theoretically 
sound, infiltration aodel that defines infiltration capacity in terms of general 
soil properties available froa soil survey maps and near-surface initial soil 
aoisture conditions derived froa soae antecedent precipitation index, a aodel 
state variable or, ideally accurate periodic reaotely sensed aeasurements. 

The approach was to use Richard's equation, describing unsaturated porous 
aedia flow, to generate a synthetic infiltration series for an array of conditions 
that could be encountered in the field within each of the SCS hydrologic soil 
groups. The results are expressed in a general fora similar to the infiltration- 
storage approach used in aany single event and continuous simulation models such 
as USDAHL (Holtan, et al., 1973), HECl (U.S. Army Corps of Engineers, 1973), 
Stanford Watershed Model (Crawford and Linsiey, 1966), and SSARR (U>S. Army Corps 
of Engineers, 1972). Infiltration-storage relationships, in which the infiltration 
capacity at any time within a rainfall event is defined as a function of the 
volume of moisture stored in some upper layer, are widely used, especially in 
continuotis simulation models, because infiltration capacity can be related to a 
period by period accounting of moisture within the SK>del components. 
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REVIE»I OF TOE INF ILTRAT I ON -STORAGE CONCEPT 
WITHIN HYDROLOGIC MODELS 


The infiltration-storage approach is examined by reviewing the functional 
form in two different Bodels: USDAHL and SSARR. The underlying assuaptions 

used in these models and all other models employing ar infiltration-storage 
approa:h are examined. 

1. EXAMPLES WITHIN TWO HYDROLOGIC MODELS 

1.1. USDAHL Model . Within USDAHL. '[.United States Department of Agriculture 
Hydrology Laboratory), an agricultural watershed model. Hoi ten's equation is 
used to express infiltration capacity ... "as an exhaustion function of available 
moisture storage in the surface horizon of a soil diminishing to a low constant 
rate of intake associated with ding strata." The equation, based on physical 
concepts as shown in Fig. 1, bv ..>rived from field infiltration data is given as: 

f - a Sj-^ ♦ fg (1) 

where: f is the infiltration capacity, in in/hr; 

is available soil porosity, undepleted by moisture, in inches; 
fc is rate of intake after prolonged wetting in in/hr; 
a is a coefficient of pore space continuity, estimated as the product 
of vegetative density of crop maturity and stage of growth. 

A graphical representation of Eq. (1) end its estimation of a typical time varying 
infiltration capacity for a soil subject to a steady rainfall is presented 
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Figure 1. Physical Concept of Holtan's Infiltration model, 
(glymph and Holtan^ 1969) 
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•s Figure 2. The aodel rests on the preaise that the sajor tone of hydrologic 
activity lies above sosie controlling depth which controls profile drainage so t.svc 
soae potential storage voluaie, S, can be conputed from total porosity and known 
or assuaed antecedent aoisture conditions. In USDAHL, primarily used in agricultural 
applications, the controlling depth, or impeding strata is considered to the 
B-horizon in agricultural soils. 

Extensive information relating the paraaeters to soil texture classes 
are available through soil surveys of the Soil Conservation Service (SCS) . 

Values of f^, the steady state liaiting drainage rate in the controlling 
B-horizon, are given in Table 1 for each SCS hydrologic soil group. The 
hydrologic capacities of texture classes are defined in Table 2. The total 
Bof storage capacity, S, is ‘he difference between moisture content 

r: i.stu«ation and wilting point. The total storage is divided into large 
pores, G, which are drained by gravity and plant available porosity, ANC, depleted 
only by plant uptake. Using this approach the moisture-retention is defined as 
three points: saturation, field capacity, and wilting point as shown in Fig. 3. 

Drainage ^nd plant uptake are determined as a function of the current storage values 
of G and ANC. The vegetative parameter, a, is defined in various landcovers in 
Table 3. An advantage of this model is that the parameters can be estimated 
easily when on SCS soil surveys are available. 

1.2 Corps of Engineers SSARR . The SSARK (Streamflow Synthesis and Reservoir Regu- 
lation) model is a hydrologic model of a river basin throughout which streamflow 
can be synthsized. The model is structured around three main coqfionents: w 

generalized watershed model, a river system model, and a reservoir regulation model. 
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FIGURE 2 Example of holtan's infiltration equation 
(unique function for particular soil) 
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TABLE l,~<41)rdro logic soil groups used by the Soil Conservation Service 


Final constant 

Hydrologic ^ infiltration 

soil group Soils included rate, (fc), in. /hr. 


A (Low runoff potential) Soils having 

high infiltration rates even when 
thoroughly wetted, consisting chiefly 
of sands or gravel that are deep and 
well to excessively drained. These 
soils have a high rate of water 
transnission. 

B Soils having aoderate infiltration 

rates when thoroughly wetted, chiefly 
aoderately deep to deep, aoderately well 
to well drained, with aoderately fine to 
aoderately coarse textures. These soils 
have a aoderate rate of water transaission. 

C Soils having slow infiltration rates when 0,05 - 0.15 

thoroughly wetted, chiefly with a. layer 
that inpedes the downward novement of 
water or of aoderately fine to fine texture 
and A slow infiltration rate. These soils 
have a slow rate of water transnission. 

(high runoff potential) 

D Soils having very slow infiltration rates 0 to 0.05 

when throroughly wetted, chiefly clay soils 
with a high swelling potential; soils with 
a high pennanent water table; soils with a 
clay pan or clay layer at or near the sur- 
face; and shallow soils over nearly inperv- 
lous naterials. These soils have a very 
slow rate of water transnission. 


0.30 -0.45 


0.15 - 0.30 
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TABLE 2. --Hydrologic capacities of vexture classes 


Texture class 

Storage 
capacity S 

Large pores 
C 

Plant-available 
porosity ANC 


(percent) 

(percent) 

(percent) 

Coarse >and 

24.4 

17.7 

6.7 

Coarse sandy loan 

24. S 

15.8 

8.7 

Sand 

32.3 

19.0 

13.3 

Loaay sand 

37.0 

26.9 

10.1 

Loany fine sand 

32.6 

27.2 

5.4 

Sandy loan 

30.9 

18.6 

12.3 

Fine sandy loan 

36.6 

23.5 

13.1 

Very fine sandy loan 

32.7 

21.0 

11.7 

Loan 

30.0 

14.4 

15.6 

Silt loan 

31.3 

11.4 

19.9 

Sandy clay loan 

25.3 

13.4 

11.9 

Clay loan 

25.7 

13.0 

12.7 

Silty clay loan 

23.3 

8.4 

14.9 

Sandy clay 

19.4 

11.6 

7.8 

Silty clay 

21.4 

9.1 

12.3 

Clay 

16. R 

7.3 

11.5 
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SOIL MOISTURE CONTENT, PERCENT BY VOLUME 
NOTE; Adoptid from 'Riprotontofivo ond Eipirimontol Botint,' oditod by 
C.Totbit ind V. Ouryvoiv, publitbid by UNESCO, 1970. 


Figure 3. Soil Moisture Constants by Texture 
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TABLE 3 

VEGETATIVE PARAMETER ■ IN INFILTRATION EQUATION f ■ 


Land Use of Cover 


Condition 

( Basal Area Rating* 

Fallow 

After Row Crop 

0.10 

After Sod 

0.30 

Row Crops 

Poor 

0.10 

Good 

0.20 

Small Grains 


0.20 


0.30 

Hay (Legumes) 


0.20 


0.40 

Hay (Sod) 


0.40 


0.60 

Pasture (Bunch Grass) 


0.20 


0.40 

Temporary Pasture (Sod) 


0.40 


0.60 

Permanent Pasture (Sod) 


0.80 


1.00 

Woods and Forests 


0.80 


1.00 


*Adjustaents needed for "weeds" and "grazing" 


(fron Glyaph and Holtan, 1S69) 
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Mlthln th« «at«rsh*d aod*}, «at«rtht4l 9rf Ifparatvd into ralativtly 

hoaoganeous units for Indapandant anaylsls bafort thay ara addad to the 

systaa. A major elaaent of tha watarshad model is the soil moisture index 

(SMI) which functions as an indicator of tha relative soil wetness and is used 

to dataraine runoff. Tha percent of tha total rainfall input that becomes available 

for runoff during some period empirically derived relationships of SMI vs. runoff 

percent (ROP). Because runoff is expressed in terms of rainfall intensity, rain* 

fall intensity may be included as a third variable in SM]<ROP relationships as 

shown in Fig. 4. The total generated runoff for period (RGP) is computed as: 

RGP - ROP X P (2) 

n 

where 1*^, precipitation for the period. 

The soil moisture index, a state variable representing current moisture 
conditions in the upper soil layer, is updated for each period by: 

SMIj • SMIj ♦ (Pjj - RGP) - ^ X KE X ETI (3) 

where SMI^ and SMI 2 are soil moisture indices at the beginning and end of 
the period; PH is period length in hours; ETI is an evapotranspiration index, 
in inches per day; and KE is a factor for evaporation reductions on rainy days. 

Nhen soil moisture conditions are approximately at the wilting point, the 
SMI is a small number which yields little or no runoff. When precipitation 
recharges soil moisture, the value of SMI increases until it reaches a maximum 
value considered to represent the soil field capacity and the runoff percent is 
assumed to approach 100 percent. 

The generated runoff for a period (RGP) represents not only rainfall 
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•xc«sf not ontoring the toil regl4e but also subsurface (upper storage) 
lateral flow and baseflow runoff. The partitioning of RCP into these three 
eoaponents is accoaplished through two additional eapirical paraactcrs: a 

baseflow index and a surface-subsurface nmoff ratio relationship. 


2. EVALUATION OF INF I LTRATICM- STORAGE APPROACH 

The USDAHL and SSARR aodels illustrate two aethodologies for estiaating run 
off for a tiae period by water budget accounting of aoisture aovenent through 
the various hydrologic conponents. In both cases, the infiltration capacity is 
described as a function of the current soil aoisture conditions (or at least a 
paraaeter related to aoisture conditions as in SSARR). This type of approach 
is advantageous for several reasons: 

1) Infiltration recovery during interaittant lulls in a rainfall event 
or between rainfall events can be accounted for through soil drainage 
and ET algorithas. 

2) Infiltration capacity can be described for non-unifom rainfall 
distributions. 

3) Coaputations are such easier than a aore rigorous swtheaatical 
approach to aiodeling infiltration. 

4) Difficult to attain soil data such as the soil aoisture retention and 
unsaturated hydraulic conductivity functions are not required to 
estiaate infiltration rates. 

jr\9 of the key assuaptions of the infiltration-storage relationship is 
that infiltration capacity can be related to soil aoisture storage without 
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r«g«rd to the v«rtlc«l distribution of watsr in ths soil profils. A rssult 
of this sssusqption is ths skistsnes of s singls-vslusd uniqus rslstion bs- 
twssn ths infiltration capacity and ths voIibm of storad Moistura in sobm 
dafinad surfaca layar within a particular soil. For axanpla, ths astistation 
of parasatars in Holtan's aquation is basad on soil prepartias or. as with 
tha ROP vs SMI raltalonship in SSARR, calibrated with rainfall-runoff data 
to obtain unique infiltration-storage function particular soil or group 
of soils within a watarshed unit. 


An axanination of tha above assunption was made by Bloomfield et. al. 

(1961) using a numerical solution of the ur.saturated flow aquation to model 

4 :onded infiltration into a semi-infinite soil column. The results of a 

series of Infiltration simulations using a range of uniform initial SMsist- 

ure conditions were expressed in terms of available storage in the soil. 

An arbitrary effective depth of soil was chosen as 200cm and the initial 

stored contents (US ) were calculated for each run asi 
o 

US^ - 200 (e. - 9 ) (4) 

o i r 

where 0^ is the uniform initial water content and 0^ is the residual water 
content following drainage. A smxlmum volume of available storage was 
defined as: 



200 (e 


0 ) 


( 5 ) 


where 0 is the saturated water content, 
max 

contents, the following relation was used: 


USt^^ - US^ ^ 41^,1 


To calculate successive store 


( 6 ) 
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where US^ it the itore contents at the end of tine period t and i* 

the infiltration voluae in tine period t«l. 

Using Eqs. (4), (S), and (6), IlcoBfield expressed the ponded infiltration 
capacity for each initial Boisture condition in teras of an infiltration- 
storage relationship. Coaparison of these relationships, shown in Fig. S, 
indicated that when the vertical distribution of soil aoisture is considered, 
the relation between infiltration capacity and soil water storage is not 
\i.<ique for the particular soil but is described by a faaily of curves dependent 
on the initial soil water storage or initial unifom aoisture content. 
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Figure 5 



Ut I I 


» ■ . ■ ^ ^ 
i *1 *4 • « • 

Uft /UftMAI I mm 

Solid lines Indicate numerical results, dashed lines 
predicted using (9) and (101. 


(BLOOMFIELD^ PILGRIM /lNr» WATSON^ 1981) 


. Nonunique aI - uS Relationships 
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CHAPTEK 3 

OEVELOPNINT AMD USE OT MCPRESENTATZVE ZNPILTAATZON-STORACE 
RELATZCMSHIPS POM MATEXSHEO NODELZNC 

1. GENERAL APPROACH 

Zf Infiltration aatiaataa aada usin9 infiItration-atora9e rataionahipa 
ara to ba ia|provad» tha work of tha Bl-joaifiald group ahowa that tha para- 
■atara of tha function ahould ba dafinad not only for particular aoil prop- 
artiaa but alao for varioua initial aioiature conditiont axiating in the aoil. 
Thia ijiprovaatant providad by incluaion of initial aoil aoiatura in a 
functional ralationahip dafinition of the iaipact of v> tarahad noiatura cond- 
itiona on pradictad rainfall axceaa volumes. 

Zt ia a difficult taak, howavcr« to define the paranctera in tarma of 
both soil properties and initial moisture conditions. For example, the 
calibration of a simple ampirically-basaO infiltration storage aquation for 
a particular watershed would require rainfall-runoff data resulting from 
an array of initial aoil moisture conditions. Extensive field tests «>*>uld 
ba required to enlarge the already existing results of inf iltrosteter data 
on field plots used to define paraateters for various soil types (such as 
Holton's equation). 

Another approach, as used in the present study, is to use the theoretical 
aquations of unsaturated moisture flow to simulate rainfall infiltration 
behavior for an array of different soils and initial moisture conditions. 
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A siailar approach waa uaad by laith (1972) to davalop a paraMtrlc In- 
filtration aquation to daacriba tha raaults of thaoratical infiltration 
aiaulationa. ly axpraaaing tha raaulta in tarwa of qanaral infiltration- 
atoraga ralationships, a aariar of rapraaantativa functiona can ba da- 
valopad. Thaaa functiona can than ba uaad to dafina tha dapandanca of 
tha infiltration-atoraga ralationahip on tha initi*' aoiatura for tha rang# 
of conditiona uaad in tha maMrical axpariaanta. 

2. DEVELOPMENT OF REPRESENTATIVE INPILTRATIOM-STORAGE FUNCTIONS 

2.1 Infiltration Modal . An array of infiltration avanta waa aianilatad 
uainq a awthanatical aiodal to daacriba rainfall infiltration into a aaaii- 
infinita aoil profila. In theory , the one-diaianaional infiltration capacity 
at tha aoil aurfaca ia baaad on tha vertical aovaawnt of aoiature in tha 
naar-aurfaca unaaturatad aoil ragiaia. Unaaturatad aoiatura flow in a non- 
awalling aoil can ba daacribad by Richard 'a aquation aai 

2 2 29 

C(9)^ - ^(9) - 1) (7) 

in which ^ la tha aoil watar praaaura haad (cai) > C(i|i) ^.a the apacific 
■oiatura capacity dafinad aa tha alopa of the aoiatura-retantion curve i 
E(9) ia the unaaturatad hydraulic conductivity function and z ia tha dapth 
balow tha aoil aurface (cb) . 

Tha aolution of Richard 'a aquation raquiraa tha apecification of boundary 
conditiona. Tha uppar boundary condition at tha aurfaca until the aoil ba- 
cooaa aaturatad and ponding occura ias 
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» - Rtf) - 1) (8) 

in which R Is ths rsinfsll flu* at th* surface. During this phase of the 
infiltration awent» the infiltration rate is equal to the rainfall rate. 

After aurface ponding occursi the upper boundary condition becoaosi 

g(O.t) - D (9) 

in v'hlch 9(0, t) is the aoil pressure head at the surface and D is the depth 
of ponded water at the surface. During this phase, the infiltration capac- 
ity, R is defined by Darcy's Law asi 

R - -R - 1) (10) 

s ps 

in which R^ is th* eaturated soil hydraulic conductivity. 

Th* lower boundary condition in th* soil column can be denoted as sosw 
depth, L, below th* soil surface where th* effect of the infiltration event 
on the initial soil conditions can be assumed negligible. Ths lo%#*r boundary 
condition becoaws 

g(L,t) ■ f(L,0) for 0< t< • (11) 

Th* initial boundary condition is ths initial volumetric stoisture content 
through the soil profile. 

An implicit finite difference formulation of Eq. (7) was solved num- 
erically subject -o the relevant boundary conditions. An implicit linear- 
isation of the non-linear functions was mad* using the predictor-corrector 
method. A complete explanation of th* formulation of th* sndel *quatlor<<, 
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th* iaplMMntation of th« nuMrlcal Bolutioiti tnd vorifioAtien of tho nuM- 
•rlcal solutions srs givsn In Nillisnin 9 

Ths thsorsticsl sodsi prsdlcts Infiltrstion rstss into s bsrs, non- 
swslliny soil bsssd on ths following input dstsi initisl soil aoisturs cond- 
itions in ths soil St th« bsginninn of ths rsinfsll svsnti s dsscription of 
ths bsslc soil hydrsulic propsrtiss givsn by ths aoisturs rstsntlon function 
snd ths unssturstod hydrsulic conductivity functloni and ths rsinfsll hysto- 
grsph. 

2.2 aisiulstsd Conditions . A ssriss of nuksrlcal sisnilstlons wars run for 
a wids vsristy of soils with distinct hydrologic chsractsrlstics. Soils %fsrs 
sslactsd thst srs rsprsssntstivs of ssch of ths fcvr SCS hydrologic soil groups 
givsn in Tabls 1. Ths SCS classification is ussd sstansivsly in siany watsr- 
shsd SKxSsls bscauss soils nan bs assignsd to hydrologic soil classss bsssd on 
tsxturs inforsuition frost routins soil survsys. As asntionsd sarlisr, ths 
basic dsscription rsquirsd for sach soil consistsd of ths soil sioisturs ra- 
tsntion curva and thv unsaturatad hydrsulic conductivity function. Ths 
raprassntstlvs soil propsrtiss wars sslactsd as followsi 

A. A soil tsxturs clsssif icstion was chossn that was rspisssntativs 
of sach SCS soil group bassd on ths dsscription givsn in Tabls 1. 
Sscauss so sksny soils fall in ths SCS B and C soil groups, two 
soil tsxturs classss wars sslactsd for sach of thsss hydrologic 
classss. Ths soil tsxturss sslactsd for sach SCS soil group ars 
givsn in Tabls 4. 
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TABLE 4. REPRESENTATIVE SOIL PARAMETERS 


( 1 ) 

Soil 


( 2 ) 

SCS 


>-■ w 
3 O 
*♦> to 

M- 

*- CO 

r» ,t 
>1 • 

m 

r» a. 
M-X 

o 

3 (O 


• r* 
r* S 
« 


^ n 
moo 
c o 
• n • 

*1 Hi 
►- Hi 
O* O H- 

^ C o 

•1 H. 


3 


m 

•d 
o 

• 9 
• 

9 

rr 


5f 


O. r» 
C 

n H 


(3) 


(4) 


CS) 


3 C 
•a H 

• 

•3 


< d 

A 4 

• o 

m n 

o O 
H- o 

3 3 

W 0 

f4 


** s 

^ 9 

o 

w o 

< 

A 

3 

a 

3 


w 

►1 

* ▼ 


(6) 

(7) 

(*) 

(9) 


Texture 

Hydrologic Soil 
Group 

f (ca/min) 
c 

K (ca/ain) 

Porosity,? ^ 

♦ jHog) 

F 

fc 

V 

Loaay 

Sand 

A 

0.013-0.019 

0.102 

0.410 

1.78 

4.38 

0.12 

0.05 

Sandy 

Loaa 

B 

0.006-0.013 

0.043 

0.435 

7.18 

4.90 

0.20 

0.09 

Loaa 

B 

0.006-0.013 

0.011 

0.45 

14.6 

3.59 

0.25 

0.13 

Sandy Clay 

Loan 

C 

0.002-0.006 

0.007 

0.42 

8.63 

7.12 

0.25 

O.IS 

Silty Clay 

Loaa 

C 

0.002-0.006 

0.004 

0.477 

14.6 

7.75 

0.32 

0.19 

Silty Clay 

D 

0. -0.002 

0.002 - 

0.492 

17.4 

10.4 

0.37 

0.26 


KJ 

KJ 


(3) Fron Table 1. 

(4) Froa Rawls et al., 1981. 

(5) (6) (7) Froa Clapp and Homberger, 1978. 
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■•cause of th« variability of hydrologic prop«rtias oven vithin 
a soil taxtur* claaa« ropraaantativa Bolatura-ratantion curvaa 
darivad from a statistical analysis of dasorptlon data (Rawls 
at al., 1961} Clapp and Hornbargar, 1978) vxe salactad for each 
soil taxtura class. An axtensiva coaipilation of hydrologic 
propartias was suanarizad and usad to astinata tha paramatars of 
tha Bioo)cs and Coray aquation by Rawls at al., (1991). Tha B' ooks 
and Coray aquation provides a reasonable astijsation of the 6 
curve for tensions lass than 50 cm (Brakensiak, at al., 1981) 
and is given as: 


X 

where is the affective saturation, is the capillary pressure 
head (cm) , is the bubbling pressure (cm) , and 1 is tha pore size 
distribution index. The parameters of an empirical power curve 

ware estimated for each soil texture by Clapp and Hornberger (1978) 
based on desorption data by Holtan et al. (1966). The power curve 
was defined as: 


♦ - ♦g(8/B,)"^ (13) 

where is the satur, ted water content or porosity; 6 is the water 
content; and both and b are empirical parameters. In addition, 
a parabolic sodification of the power curve was proposed to account 
for gradual air entry near saturation. The functions resulting from 
Bqs. (12) and (13) were compared for specific soil texture 
classes and found to be vry similar. The power curve, given by 
Eq. (13) was used with representative parameter values in 
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Table 4 froa Clapp and Hornbar 9 cr (1978) . The resulting aolature- 
retentlon curves are given in rig. 6. 

C. Representative unsaturated hydraulic conductivity functions were 
described by a swthod proposed by Casipbell (1974) . The unsaturated 
hydraulic conductivity function is described based on a linowledga 
of the power curve describing the sioisture retention function and 
and estiawte of K^. The equation is given asi 

*(♦) - K^(e/e^)^‘^^ (14) 

where K(i|i) is the hydraulic conductivity at a specified soil pres- 
sure, ill and b and are defined as in the power curve ^iven in 
Bq. (13) . An estiMte of a representative value of for each 
soil texture is difficult because of the large variations of this 
property even within a soil texture class. A set of aean values of 
saturated hydraulic conductivity values for each soil texture class 
given by Rawls at al. (1981) were used in this study. 

For each soil given in Table 4, infiltration sisiulations were run for 
a range of initial soisture conditions that %«ould be expected to occur under 
field conditions. Soil soisture conditions that were tested ranged from near 
the wilting point to near saturation. Moisture conditions at wilting point, 
field capacity and saturation, defined lay the snisture retentio.'! curve are 
given in Table 4. A sussnary of the siaulated initial conditions is given in 
Table 5. The vertical distributions in each initial Moisture profile is as- 
sixsed to be unifora. The ixqalications of this assuMption are discussed in a 
later section. 
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FIGURE 6. REPRESENTATIVE MOISTURE-RETENTION CURVES 
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INITIAL NOIITUU CONDITIONS fINUUTEO 
FOB BACH MPBESENTATIVB 901 L 


Soil 

SCS Group 

Initial Volusatric Noittur# 

% Saturation 

Lossy 

A 

0.24 

59 

Sand 






0.38 

93 


Sandy 

B 

0.28 

64 

Loaa 


0.40 

92 


Loan 

B 

0.19 

40 



0.22 

49 



0.28 

62 



0.32 

71 



0.36 

80 



0.40 

89 

Sandy Clay 

C 

0.22 

52 

Loan 






0.30 

71 



0.38 

90 

Silty Clay 




Loan 

C 

0.25 

52 



0.35 

69 



0.45 

94 

Silty 




Clay 

D 

0.34 

69 



0.40 

81 



0.46 

93 


26 




Of POOR 

2.3 D»vloB— nt of Infittr>tion-»tor«Q« Functions . The results of each 
nuMrical alaulation dafinad infiltration capacity aa a function of tiaw for 
a particular aoil, initial aoiatura condition and rainfall avant. An ax- 
aapla of rasulta of thaaa aimulationa is shown in 1*19. 7. whara tha infiltra- 
tion capacity of a loaa soil with an uniform initial aoiatura content of 0.26 cm/cm 
is shown during thraa rainfall intansities. 

Rasulta such as thosa shown in Fig. 7 ware axprassad in an inf iltration- 
storaga ralationship as follows: 

A. An incraawntad tiSM period was defined for period by period 
accounting of soisture axivancnt. This tiae period should corres- 
pond to the time period used to describe the hyetograph in the 
particular siaulation model in which the infiltration storage 
relationship would be used. A time period of 15 siinutes was chosen 
as a representative value to be used in this study. However, in- 
filtration-storage relationships could be developed from the simula- 
tion results for any time period of interest. Tl.e time period is 
denoted as At. 

B. Two quantities were defined as follows. The depleted storage, 

AS^, was defined as: 

AS^ - ri^ (15) ■ 

whara AS^ is the total voluaatric change in moisture storage 
from initial s»isture conditions at time t} and 11^ is the total 
voluBW of infiltration frcai the rainfall avant at time t. The in- 
filtration volume ovar a particular tima period, Ai^, is 


Q’jALrTV 
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INFILTRATION RATE, I (cm/MIN) 
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REPRESENTATIVE LOAM SOIL 

- Kj^O.Ol cm/min ■ 0,25 in/hr 
— POROSITY ■ 0.45 cm/cm 



FIGURE 7. TIME VARYING INFILTRATION CAPACITY 
RESULTING FROM NUMERICAL SIMULATIONS 
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d«fio*d aai 




(16) 


wnara &S.«nd arc tha d«platad atora 9 # at tha baginnin^ and 

V oX, 

and of tha tiaa pariod» raapactlvaly. 

C. Uaing aquatlona (15) and (16), tha Infiltration capacity ovar a tiM 
atap, ^ axpraasad as a function of t)ta daplatad ator- 

aga at tha ba 9 innin 9 of tha tiac stapAS^. An axanpla of this 


transfonsation is shown in Table 6 for tha rasnlts of a particular 
siaulation shown in Fig. 7. From tha numarical rasults, the cum- 
ulative Infiltration volume at 15 sdnuta time incraaants are datar- 


sdnad. Using Equation (15), can ba equated to the cumulative infil- 
tration volumes as shown in Table 6, cols. (1) and (2). Using Eqjstion 
(16), infiltration voIusms ovar the time step are defined (col. 3). 


Tha infiltration-storage relationships resulting from tha rainfall 
simulations for a loam soil with an initial moisture conditions of. 0.36 c 3 n/cm 
(Fig. 7) are shown in Fig. 6. Whan the H ~ tS relationship is used to 
describe tha inf i Itration capacity, a unique function, auch as shown in Fig. 6, 
describes the infiltration behavior of a soil with a specific initial sioisture 
condition ragardless of the rainfall rate. Of course, prior to surface 
saturation and ponding, the infiltration rate into the soil would k>e limited 
by the rainfall rate as sho%m in Fig. 8. 

The fact that the Bl - relationship for a specific soil and initial 
snisture condition is independent of the rainfall rate is significant because 
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TABLE 6 

INFILTRATION- STORAGE REUTIONSHIP 
DERIVED FROM RESULTS OF A NUMERICAL SIMULATION 


(3) 

I (c») 


( 1 ) 

t(«in) 


( 2 ) 

S^(c«) 


0 

0 

IS.O 

0.786 

30.0 

1.325 

45.0 

1.730 

60.0 

2.077 

75.0 

2.387 

90.0 

2.671 


.786 

0.S39 

0.405 

0.347 

0.310 

0.284 
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INFILTRATION VCLUME^ Al (cm/15 MINUTE PERIOD) 
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FIGURE 8 INFILTRATION-STORAGE RELATIONSHIPS FOR A 
LOAM SOIL (6j ■ 0-28) 
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it indicates that this function, dsrlvsd froai nuawrical aijRulationa for a 
apscific rainfall rats, can bs used to satiMts infiltration capacity during 
any apscifisd rainfall avsnt. Ths hypothesis fonss the basis for a sinplc 
proesdurs to bs presented in the next section. 

The tl - AS functions resulting froai a range of initial eioistura conditions 
in the loaai soil are shown in Fig. 9. The family of curves are similar to 
those presented by Blooaifield et al. (1981) in Fig. 5 showing the isipact of 
the initial aoistura condition on the structure of the infiltration-storage 
function for a particular soil. The functions presshted in this study are 
different in structure from those developed by Bloomfield however, because 
of the way that soil water storage is defined. As outlined earlier, Bloom- 
field used equations (4), (5) and (6) to express I as a function of used 
storage, US. The quantity, US, is a function of both initial moisture cond- 
itions and a specified depth of evailable storage and, thus, the curves are 
unique for the depth of storage that was assumed. 

The curves shown in Fig. 9 are defined in a more general form by ex- 
pressing infiltration capacity, £1, in terms of depleted storage which is, 
by definition, the cxssulative infiltration volume. In this format, the 
curves can be applied to any specific watershed where the depth of avail- 
able etorage in the upper soil Isyer stay vary. 

The dl - dS functions are given in Appendix A for each of the six soils. 
Functions ar. given for time periods of 15 adnutes and 5 minutes. 
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Fr-URE 9. INFILTRATION CAPACITY AS A FUNCTION OF DEPLETED SOI.. MOISTURE STORAGE 



0 . 1 ** • » » 1 » ~ » "» ~ 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 ^.0 


DEPLETED STORAGE FROM INITIAL CONDITIONS «s(cm) 

(cumulative infiltration) 
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3. use Of ZMPXLTRATION - STORAGE FUNCTIONS 
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3.1 Exmi\p1« problwn . The approach xa bast lllustratad through tha uaa of 
an axanpla givan in Tabla 7. Tha illuatratad rainfall hyatograph occurs 
on an SCS B hydrologic aoil with an initial voltvatric aoiatura cond- 
ition of 0.28 cm/cm. Tharafora, curve 4 for a loam aoil (Fig. 9) is uaad 
to aatiaata infiltration capacity throughout the event. Tha procedure is 
outlined as follows! 

1) Coluatns (1) and (2) are defined by the rainfall hyatograph for 
tha specific problan. Column (1) is the time interval by 15 mir- 
uta increments. Colav) (2) is the volume of rainfall occurring 
during each 15 minute incrastent (for example, 0.26 in/hr x 15/60 x 
2.54 as/in ■ 0.165 cm). 

2) AS, in Column (3), is tha volume of available storage in tha soil 
colusin that has been depleted at the beginning of each time increment. 
By definition, AS > 0 at the beginning of the rainfall event be- 
cause no storage has been depleted. 

3) AI, in Column (4) represents the total volume of rainfall infiltration 
during the time increaient. Al defined as a function of AS (Col. (3)) 
by using curve 4 in Fig. 9. 

4) Q /n Column (5), is the rainfall excess volume over tha time 
increment (A R - 61) and ZQ (Col. (6) is the cumulative rainfall 
excess volume. 

5) At the beginning of the 0.15 minute increment, no storage has been 
depleted. Using curve 4, the infiltration capacity for 6S ■ 0 
exceeds 0.6 cm, so the entire voIusm of rainfall, 0.165 cm, in- 
filtrates into the soil. 
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TABLE 7 

EXAMPLE PROBLEM: USE OF Al - AS CURVES 


E»ijpl« i Estinate Infiltration for Following Conditions: 


1. Stora hyatogrsph 



2. Soil properties: SCS soil group B 

(assuae Loaa soil) 

S. Initial aoisture conditions: 0.28 ca/ca 

4. Assuae unliaited depth of storage 

Procedure: Use representative infiltration-storage functio.i to describe 

infiltration capacity for this soil and initial aoisture condition. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Tiae (ain) 
Interval 

AR/15 Bin 
(ca) 

aS(cb) 

Al/lS Bin 
(ca) 

Q(cm) 

tQ(ca) 

rQ(cm) 

numerical 


0-15 

0.165 

0.0 

0.165 

0.000 

0.000 

15-30 

0.165 

0.165 

0.165 

0.000 

0.000 

30-45 

0.240 

0.330 

0.240 

0.000 

0.000 

4S-60 

0.240 

0.570 

0.240 

0.000 

0.000 

60-75 

0.630 

0.810 

0.555 

0.07S 

0.07S 

75-90 

0.630 

1.365 

0.400 

0.230 

0.305 

90-105 

0.480 

1.765 

0.342 

0.138 

0.443 

105-120 

0.480 

2.107 

0.310 

0.170 

0.613 

120-135 

0.315 

2.417 

0.285 

0.038 

0.643 

135-150 

0.315 

2.702 

0.270 

0.040 

0.688 

150-165 

0.240 

2.972 

0.240 

0.000 

0.688 

165-lSO 

0.240 

3.222 

0.240 

0.000 

0.688 



RAINFALL VOL-4. 14ca-l .63 in 


RAINFALL EXCESS-0. 69ca 







Of POOR QUAir^ 


6) At th« b* 9 lnnin 9 of th« lS-30 mlnutf; ^t^.'lod, 0.16S c» of ttorag* 
has boon doplotod • rooult of tho infiltration dviring tho 0-lS 
adnuto poriod. Again curvo 4 indicataa that for AS • 0.165, 

tho infiltration capacity oxcooda tho rainfall rata, and tho ontiro 
rainfall voluato infiltratoa. 

7) Tho abovo procoduro ia ropoatod for oach tiao incraaont by adding 
tho infiltration volvato, A I, froa tho prcvioua tiao atop to the 
cuBulativo doplotod atorago. As. The AI - AS curve ia then uaod to 
define the infiltration capacity for the current tiaj atop aa a 
function of As. 

8) In thia exaatple, all precipitation antera the aoil until the 60-75 
alnuto increment. At the beginning of thia increment, AS > 0.910cm 
and the reaulting Al _,^ (from curve 4) ia 0.555 cm. Therefore, the 
infiltration rate ia now limited by the aoil capacity and the rain* 
fall exceaa (Col. (5)) ia 0.630 - 0.555 - 0,075 cm. 

3.2 Coopariaon with Mumerical Simulation . The reaulta uaing thia procedure 
were compared with the infiltration capacity reaulting from a numerical aol- 
ution of Eq. 7 uaing the apeclfic rainfall boundary conditiona of thia prob- 
lem. The cumulative rainfall exceaa volumes at the end of each time increatent 
are given in Column (7) of Table 7. A compariaon of Columns (6) and (7) Indicates 
that the use of the AI - AS rclatlonahlp pro"’ded results in very close agree- 
ment with the numerical simulation of a specific rainfall event. The infil- 
tration rate cruves from the numerical simulation and the inf ilvration- 
atorage method are shown in Fig. 10. 


INTENSITY (in/hr) 


ORIGINAL PAOZ 
Of POOR QLIALrrY 

I 

TOTAL RAINFALL VOLUME “ 1.630 IN. 

PREDICTION OF RAINFALL EXCESS 

- INFILTRATION STORAGE METHOD Q “ 0.271 IN. (16.6Z) 

- NUMERICAL SIMULATION Q - 0.269 IN. (16.5X) 



Figure 10. Comparison of Results Using Infiltration-Storage 
Hethod and Numerical Simulation 
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4. oowiiomATioii iM Th£ on p yiiriLTiuTiqn -|ypyi^i cimvyy 

Th« Al - A6 curves were developed froei nunerlcel elauletione for e 
eljc>lified eoll eyetem. In order to epply these curves in practical water- 
ahed Bodeling efforts, it Is necessary to consider how to adjust the curves 
to aiaulata conditions where the assumptions used in developing the curves 
■ay not be valid. Such conditions might include the followlngi a) some 
specified raxisnia depth of available soil storagei b) watershed soieture 
accounting SKxlals with varying simulation time periodsi c) a non-uniform 
vertical distribution of initial eoil moisture in the soil column* and d) 
variability of important soil hydraulic properties such as the saturated 
conductivity and porosity from the assusMd representative values f> a partic- 
ular soil texture. 

4.1 Maximum Depth of Storage . In most continuous simulation models, the 
infiltration capacity is related to some defined mixisvan available storage 
reflecting both the initial moisture deficit and a finite depth of soil 
as the estimated depth above aome impeding soil layer of lower hydraulic 
conductivity (the A horison in agriculture). In some models, the maximum 
depth is implicitly estimated in the form of an upper soil storage parameter 
that is defined through empirical approach. 

As mentioned earlier, the tl - AS curves were developed from nusierical 
simulations of infiltration into aoil columns of unlisdted depth of available 
storage. Using this approach, curves developed have a s»re general nature 
instead of being structured for a specific depth of storage. However, the 
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curvci could b« applied to • watershed condition haeing a specified depth 
of available storage by defining the SMixiausi available storage asi 

^®«ax " ^®s ■ *1^** 

in which d is ths dspth of available storage to surface infiltration. The 

infiltration capacity of the soil is then defined by the - tS function 

until ds, the depleted storege, approaches at this point, the svail- 

abls storage in the upper soil sons is exhausted and infiltration capacity 

Bust be defined in sons other Banner consistent with the structure of the 

overall Bolsture accounting Bodel. In Holtan's equation, for instance, the 

infiltration capacity would be f , the steady state liaiting infiltration rata, 

c 

whan surface storage becoaas exhausted. Zn SSAMt, when SHI reaches a BaxiauB 
corresponding to fisld capacity, all rainfall becoaes either surface or sub* 
surfacs runoff. From a physical standpoint, once the surface soil storage 
is exhausted, the infiltration rate would be governed by the saturated con- 
ductivity of the surface soil layer of the iapeding layer controlling drainage 
froB the surface soil storage layer. 

Figure 11 shows how the general AI - AS curves for a loan soil could be 
Bodified to represent a field condition where the available storage to a 
surface layer of seven inch thick soil overlaying an essentially inpervious 
lanse. Hhen the storage of the seven inch layer of soil is exhausted, the 
infiltration capacity goes to sero or becoaes controlled by the Bore im- 
peding soil. Using this type of axidification, allows for infiltration cap- 
acity tit be described in tens of soil texture and the depth of available 
storage which considers basin geoenrphology and physiography. 
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4.3 Sijmlation Tl— P«rlod . Th« & I - IB ralatlonahlp is defined in terms 

of • specific tire incresMnt. In sumy instences however, it mey be desirsble 

to estiaate infiltration capacities based on sos« other accounting period 

used within a particular watershed siodel. As Bloomfield, et al., (1981) pointed out, 

it would be incorrect to merely use the relationship from one accounting 

period to estimate infiltration vol\ates for another accounting period because 

the tl ~ AS relationship is unique for the specific accounting period. The 

isf>act of the time period is Illustrated in Fig. 12 where the infiltration 

capacity is expressed cm/hr for 15 sdnut-. and 2 minute increments. 

The AI - AS functions were expressed for two time periods of 5 and 15 
minutes which were chosen tr be representative of the time intervals in 
simulation models. In addition, time intervals of greater duration could be 
handled by breaking each increment into one of the t«K> increments for which 
Al - As function was defined. For example, if the simulation model required 
hourly precipitation and evaporation data and rerformed hburly accounting of 
t>oisture SK>vement, the rainfall hyetograph given in hourly increments could 
be divided into 15 minute increstents and the appropriate function could be 
used to estimate infiltration for successive 15 adnute periods throughout the 
rainfall event. Even if the overall simulation sodel uses time movements 
of one hour or greater, it would be advantageous to estimate Infiltration 
capacity with norm accuracy over shorter-time intervals. 

4.3 Non-uniform Initial Moisture Profile . As mentioned earlier, the 

4i - AS relationships were developed for initial soisture conditions vertically 

uniform throughout the soil registe. Without using such an assumption, an 


41 







ORlGiNAL PAG? !3 
OF POOR QUALfTY 


infinite nunber of initial conditions could have been ueed. Thus, the 
AX “ AS curves define t.'ie infiltration capacity for a particular initial 
•oieture condition throughout the entire soil coluain. Moisture conditions 
are not necessarily unifora in the colunui and two questions arise. First, 
just how isqmrtant are reaotely sensed BeasureBents of soil soisture in the 
top few centiaeters of the soil given the uncertainty of the SK.isture con- 
ditions deeper in the profile which cannot easily be Beasured? Secondly, even 
if such inforaation on the distribution of initial soil Boisture conditions 
existed, what value would it have in a practical sodeling approach? 

In order to answer the above questions we Bust consider the sensitivity 
of the infiltration procees to initial siolsture conditions at various depths 
in the soil colusm. The probleai is complex because of the variety of soil 
properties, initial sioisture conditions and rainfall conditions. A rational 
exaalnation of the theoretical infiltration process Indicates the following: 

1) The infiltration process is Bost sensitive to initial sioisture 

conditions near the surface because: a) the time to surface pond- 

ing and initiation of runoff is controlled by the tijne required to 
aaturate storage near the surface; and b) the infiltration capacity 
during ponded infiltration is controlled by the Boisture gradient 
at the surface (Eq. I'O . 

2) The ispact of the initial Boisture conditions on the infiltration 
process decreases with depth because the impact on surface ponding 
and Boisture gradient decreases. 

3) Obviously, the deeper the wetting front Boves into the soil, the 
greater the impact of Initial Boisture conditions at soeie depth 
below the surface. 
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Th« abov* obsarvationa wara axaainad fox a particular aoil ualnq 
nuMrical alanilatlona. Four initial aoiatura profilaa all with a aurfaca 
soil aoiatura of 0.28 cm/cm ara shown in Fig. 13 at profila 1 raprasants 
tha caoa of a unifora initial aoiatura condition; profila 4 raprasants a 
aora rapid variation of initial aoiatura condition naar tha surfacai and 
profilaa 2 and 3 raprasant intaraadiata caaaa. Infiltration capacity was 
sianilatad for aach initial aoiatura profile froei a rainfall of 0.75 in/hr. 

Tha rasultlng rainfall axcass volumes for aach initial aoiatura profila ara 
shown in Fig. 13 b at various stages of the rainfall. The results indicate 
thara is no significant diffaranca batwaan rainfall axcass volumaj rasulting 
from profiles 1 and 2 for any of the rainfall durations. The differences 
resulting from profile 1 and both profiles 3 and 4 were very snail for the 
short duration events but incraasad as the rainfall duration increased and 
the wetting front moved further into the soil. This is expected because 
the initial moisture content will only affect tha infiltration process near 
the surface to a depth which has been reached by the wetting front. The 
longer the infiltration event progresses, the greater the depth of initial 
moisture will impact the process. Still, the rainfall excess produced by 
profile 4 for the three hour event was less than 10% greater than that 
estimate with profile 1. 

This example rapreaants one specific case and in no way should be 
used to avsJ.uate the iag>act of initial moisture at various depths for the 
ganar'&l case of all soils and initial conditions. The example discussed for 
a loam soil, representative of an SCS B soil, indicates that surface measure- 
ments are important even though tnere is uncertainty concerning soisture 
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conditions at greater depths. Obviously, additional Information on conditions 
throughout the profile could lsg.rove Infiltration astiawtes but these lai- 
provesMnts may not be elgnlflcant given other uncertainties In natural water*- 
aheds. 

kn eaaadnatlon of the physical principles driving the Infiltration 
process indicates that for a variety of soil textures, a knowledge of the 
initial Bolsture conditions below a near aurface region of the top 2-S cm 
stay not be necessary for a good estimate of Infiltration and resulting 
runoff. In aolls of high permeability, such as sand and sandy loam (5CS 
A or B soils), the moisture front SMy advance Into the soil rapidly and be 
affected by initial moisture conditions at greater depths In the soil. How- 
ever, in these soils, the hydraulic conductivity Is large and dostinates the 
process resulting In a small runoff volume regardless of the initial SK>isture 
conditions. In soils with lower permeability, clay and silt loam, for ex- 
ample iSCS C and D solid), the moisture front seves very slowly Into the 
soil because of the low hydraulic conductivity. Because the moisture front 
advances slowly. Initial solsture conditions at greater depths do not affect 
the process. In thsse soils, the runoff potential Is high and greatly Im- 
pacted by the initial moisture conditions near the surface. 

In some cases, data stay be available describing the distribution of 
soil moisture with depth and the question arises as to the value of such 
data in practical modeling. Such lnfo:nnatlon as vertical distribution may 
be available from: 1) store sophisticated remote sensing tectiniques with 
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longer wevelengthsi 2) •stenslv* aurfact aaaauraaMnts to graatar depths 
ualng theoretical aodals driven by soil texture proparitea and clieatic 
datai and 3) extensive fiald BsasuraaMnts. 

Tha additional data could ba incorporatad into tha inf iltration- 
storaga approach if dasirad. Tha tl - tS relationship resulting froei 
profile 4, (Pig. 13) is shown in Pig. 14. In tha early stages of infil- 
tration, tha dX - AS curve follows tha curve corresponding to ■ 0.28 
and as As increases the curve approaches tha curve corresponding to • 0.32. 
Me Bight attasipt to represent this unique curve as a function of the curves 
representing the t««o uniforv Boisture conditions. One sixipla approach 
is outlined in Pig. 15. Starting with an assusMd initial soisture condition 
(Pig. 15 a), the profile is represented as a series of superimposed rectangular 
volusies as shown >n Pig. 15 b. The sMxixiuB storage is calculated for each 
volueie as the Boisture deficit (0^ -8^) x the depth of the volume eleaent. 

Tha A1 - As curve corresponding to the initial soisture content of the 
first voluBe is used to calculate infiltration capacity until the BaxiBum 
storage has exceeded the first volume. The AX - AS curve for the initial 
Boisture condition of the next volume is then used to estimate capacity 
until the Baximvxn volume of storage has been reached. Tha A X > A S curve 
is first artered at a AS corresponding to the voluaw already depleted by 
the preceding volimie element (the area of overlap) . Tha process could con- 
tinue for as many volume elements as required. Tha resulting approximation 
of the AX - AS function would be sixdlar to that shown in Pig. 15 c. The 
method assimws that one volusw alaaMnt is coapletely filled before soisture 
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Voluaerric Soil Hoisturt OR’f*»NAL PAwE tS 



(15«) Assuaed initial aoisture 
distribution 


(l^b) Siaplifiod Reprasentat io 
of Storage Volume 



(15c) a 1 - Function 
Approximated for Non-Uniform Profile 


FIGURE 15 . USE OF AI -AS CURVES TO REPRESENT NON-UNIFORM INITIAL CONDITIONS 
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K>v«i into th« n«xt olOMnt. This •pproxlMtion Mould b« momt rcasonabl* in 
•oils in which tho ■oisturo wotting front is vsry stssp« such ss sand and 
ioasi, as opposad to silt loasi or clay <oat snistura-rstsntion curvss in 
riy. 6). In addition, this approxiaation would ba scat valid in a dry initial 
condition whara tha watting front tands to be staapar than in a wat initial 
condition. 

An axaapla of tha abova procadure is shown in fig. 16 for a particular 
non-unifoni initial aoiatura profile. The runoff volwe resulting frost a 3 
hour rainfall of 0.75 in/hr (2.25 in. total rainfall) was calculated as 0.77 
inches. This runoff voluae corresponds closely to tha vnlu»e of 0.76 inches 
resulting froa a nusiarical siaulation of the actual initial profile. The 
runoff voluues resulting froB only using the 63 ~ bS curves for uniform initial 
Boisture profiles of 0.28 and 0.32 CB/cm were 0.66 inches (11.6% error) and 
0.63 inchas (7.6% error), respectively. Although the axaxple indicates that 
the curve adjustaents will wor)c, the cifferences caused by the profiles in- 
vastigated ware relatively SBall. Because of the materous other uncertainties 
in a natural watershed, it is doubtful that tha lisdtad gains accaeg>lished 
by tha more coBplax adjustaants are worth the effort. 

4.4 Variation of Soil Hydraulic Properties . Tha use of the I - S curves 

assiates that the particular hydraulic properties of the soil (K , 6 , K - «|», 6 - i(t) 

s s 

are siailar to the representative values used in this study to derive the curves. 
Unfortunately, field studies and laboratory tests have indicated that important 
hydraulic properties such as K can vary significantly within soils of tha same 

O 

texture class. Field data presented by Nielson at al., (1973) indicated that 
the coefficient of variation (CV) in could ba 50% greater within fields of 
generally hoBogeneous soil texture. The effect of variation in on the 
AI - AS relationship for loam and silty clay loam soils is shown 
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Problea: 


EitiMte infiltration voluue 

• SCS B soil 

• non-unifom initial condition 
’ rainfall 
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in ri 9 B. 17 and 18, raapactlvaly. For ««ch soil, K, was variad froa tha 

• 0 

raprasantativa valua by • ona atandard daviation aaauning a noraal dlstribu* 
tlon with CV of 50%. Thus, in using tha raprasantativa dl - dS curvas to 
rapraaant a particular soil within a wstarshad, wa might axpact significant 
arrors in tha astiMtion of infiltration capacity dua to tha uncar tainty 
of thasa hydraulic propartias. 

Tha variation of soil hydraulic propartias within a soil taxtura 
class must alwsys ba racognixad whan attasipting to pradict runoff potantial 
of an ungagad watarshad basad only on SCS soil groups or soil taxtura classas. 
Howavar, in SMny instancas tha raprasantativa &1 - tS curvas smy prova 
usaful daspita tha problam of variability of soil hydraulic propartias 
within soil made for larga ungagad watarshads whara soils ara catagorlzad 
according to SCS soil groups and surfaca moisture conditions scnitorad by 
ramota aansing. Nhan a watershed is gaged with historical rainfall- 
runoff data, or when more detailed information on soil hydraulic properties 
ara available, a different method smy ba more feasible. 

A possible approach to incorporating tha variability of hydraulic 
paramatars t«ould ba to generate a family of dX - tS curvas using a Honte 
Carlo sisnilation. Using this approach, the dl - dS relation for a part- 
icular soil taxtura having specific initial snistura conditions would ba 
axpressad as a probability function rather than ona specific curve. Ideally, 
such a probability function «K>uld Include the effects of variation of hy- 
draulic properties within the soil taxtura class. The number of simulations 
required to develop such a probability function for each soil texture and 
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initial aoiatura condition would bt axtroaely larfe. In addition, it would 
ba vary difficult to rapraaant tha cumilatlve affects of variation of nan/ dif- 
ferent hydraulic properties. 

A Bore practical approach to including the effect of varying soil properties 
within each texture class sight be to use the - tS curves in a form siailar 
to that found ill the Stanford Watershed Model (Crawford and Linsley, 1966). In 
SNH IV, areal variations of infiltration capacities are defined by plotting 
a cuaulative frequency distribution of infiltration capacity as shown in Fig. 19. 
In theory, this curve would result if a large nuaber of siaultaneous infiltroaeter 
aeasureaents were aade and plotted to show the percentage of watershed area with 
an infiltration capacity equal to or less than the aeasured values. Within a 
field of hOBOgeneous soil textures, the variability could be largely attributed 
to variations in soil hydraulic properties. If one wished to consider the var- 
iability of K occurring within a texture class when estiaating the infiltration 
capacity using the AI - AS functions, one approach could involve the following 
steps: (shown graphically in Fig. 20) 

1) Froa field data, estiiute a theoretical distribution of K values 

s 

within a particular soil texture class (Fig. 20 a). 

2) In addition to the A I - AS function already generated using the 
representational value, generate at least two aore Al - AS functions 
with varied from the representative value (Fig. 20 b) . 

3} For a range of AS values, Al is plotted for each K value as a 

function of the fraction of the watershed unit assxaed to have that 
particular value of Kg (Fig. 20 c) . 
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Pisisture 

Supply 



Equal to or less than the Indicated Value. 


(from Crawford and Linsley, 1966) 


FIGURE 19 . INFILTRATION APPROACH IN STANFORD WATERSHED MODEL 


56 



ORlGIf 
Of PC 

(a) AtauB* dittributlon of 
within aoil taxtura. 


(b) Generate Al - AS curves 
for variation of K^. 


(c) For a spacific AS , axpress 

variation of Al over watershed 
unit (following approach of 
Sm IV). 



FIGURE 20. CONSIDERATION OF SPATIAL Vi 
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4) Th« ruinftll aoiiturt supply AR is partitioned into infiltration voluae 
(the area under both the &R and dl curve) or runoff voliaae (the resiain* 
ing area under the aR curve) fflg. 20 d). 
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CHAPTER 4 

APPLICATION OF INFILTAATION-STORACE METHOD IN WATERSHED 

Tht infiltration«stor«ge aathod presented in this study could be used 
to laproye estivates of the soil infiltration capacities which are crucial for 
the accounting of voisture voveaent throughout a co^>rehensive watershed model. 
A generalized flowcliart of hydrologic processes occurring within a watershed, 
in Fig. 21, shows how the infiltration component affects the flow of moisture 
through the system. Regardless of the particular model structure, the division 
of precipitation into surface and subsurface processed is an important element 
in estimating surface runoff and updating soil moisture conditions that will 
affect future events. 

The infiltration-storage method may be useful within a new generation of 
watershed moisture accounting models structured specifically to maximize the 
benefits of remotely sensed data. This type of model is one ultimate goal 
of current research efforts at RSSL and might be configured as shown in 
Fig. 22. The key components of the model would include: 

1) A geo-referenced grid cell data base containiny spatially distributed 
soil texture classifications obtained from soil survey maps; land- 
cover information from Landsat and aerial photography; and other 
relevant spatial data. Moisture accounting could occur at a spatial 
resolution equal to that of the data base. 

2) A continuous moisture profile accounting model for some defined upper 
soil store zone affected by surface processes. The mo<tpl component 
would predict vertical soil moisture distributions subject to ET and 
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FIGURE 22 . REMOTE SENSING-MODELING LINKAGES FOR COftTINUOUS 
MOISTURE ACCOUNTING AND INFILTRATION ESTIMATES 
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d««p percolation and would b« driven by periodic reaotely aenaed aur* 
face toil t«ijperaturea and •oieture aeasureoMnis and other available 
cliaatlc data. 

3) The rainfall infiltration coaponent which would provide the crucial 
partitioning of rainfall into surface runoff and toil aoisture re* 
plenishiMnt. The aodcl component would predict tine varying infiltration 
capacity given a tlae varying rainfall hyetograph, a soil texture class- 
ification, and initial soil aoisture conditions. Initial soil aoisture 
conditions would be represented by either an Initial toil aoisture pro- 
file accounting aodel in which the initial surface aoisture is updated 
by periodic reaotely sensed aeasureaents. 

The - &S curves could be used to provide the infiltration coaponent. The 
approach is theoretically sound yet coaputational ly staple enough to calculate 
infiltration rates for aany spatial units representing various soil texture 
and Initial aoisture conditions in the watershed. 

It is iaportant to consider how existing hydrologic aodels could be 
Bodified to gain the full benefit of reaote sensing technology. These pos- 
sible aodifications are iaportant because a significant effort has already been 
invested in the developaent of accepted aodels and the aissions of aany organiz- 
ations revolve around their use. In addition, the acceptance and Justification 
of a new data source such as that offered by reaote sensing will aost likely 
coae as the utility of such data is deaonstrated in current aodrling require- 
aents. Incorporating the AI - AS curves into current aodels would provide a 
siaple, yet theoretically sound infiltration aodel that will logically incon>orate 
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th« lapact of ontocodont soil aoisturt conditions on vatorshed rosponst. The 
operation of this tjrpe of watershed «odel would interface perfectly with the 
anticipated capabilities of remote sensing aeasurewents of surface soil smisture. 


The effort involved in implementing the approach within existing models 
would vary considerably depending on the particular model structure, for 
example, within USDAHL and other simulations models using Holtan's infiltration 
equation, it would be relatively easy to substitute the representative 
tl > bS curves in place of the infiltration function described by Holtan's 
equation. The reason is that both Holtan's equation and the representative 
tl - tS curves relate infiltration capacity to the sane physical quantity, the 
volumetric soil moisture storage in a defined upper soil zone. The tl - tS 
approach could be impler^Ded using the FORTRAN subroutine given in Appendix B, 
which calculates the infiltration capacity using tabular values of the 
representative tl • tS curves given in Appendix A. Another approach would 
be to fit the function described by Holtan's equation to the representative 
tl ~ tS curves. The parameters of Holtan's equation would then be a function 
or soil properties and initial moisture conditions. Bloomfield et al., (1961) 
demonstrated this process by fitting two other infiltration storagei equations 
to theoretical infiltration-capacity curves for a specific soil. 


The use of the representative dl - AS curves within the SSARR model 
would be more difficult because of the scdel structure. The moisture accounting 
position of the model, shown within the box in Fig. 23, consists of three 
tabular functions vhich are used to separate infiltration and suS^quent 
percolation and Interflow from immediate surface runoff. The runoff percent 
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(ROP) vt soil Moisture index (SNl) table defines the eaouni * rainfall 
ifhich Infiltrates to tatisfjr soil Moisture deficits. Nhen the soil Moisture 
content reaches field capacity, the Moisture deficits are considered replenished 
and SMI reaches a MaxiMua value. The reautining rainfall, defined as runoff 
generated period (RGP}» is than subjected to a baseflow loss considered to 
represent percolation froM the upper soil zone to a lower storage zone. The 
loss is based on a table relating percent of surface runoff becoming baseflow 
(BFP) to a baseflow Infiltration Index (BII). The remaining generated 
surface runoff (RGS) is further divided ir.tu direct surface runoff (R^> and 
infiltrated interflow or subsurface runoff (RSS) according to a table 
describing the surface- subsurface separation as a function of the total 
runoff, RGS. The strategy is for these tables to be adjusted by the user to 
fit the iiodel response to rainfall-runoff historical data for a particular 
watersheil. 

The structure of the Moisture accounting component within SSARR does 
not permit the direct use of the bl - AS curves. For example, the A] - &S 
curves could not sisqply be substituted for the ROP vs SMI table because the 
ROP vs SMI table only defines the rainfall infiltration required to replenish 
soil moisture deficits below field capacity. Additional surface infiltration, 
becoming deep percolation and interflow, is considered in other tables. 

The Al - AS curves, on the other hand, predict the total infiltration 
capacity until the surface soil Moisture zone is entirely saturated. The 
infiltrated moisture would be subject to soil moisture replenisJiffent, deep 
percolation and interflow. 
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Furt^^r, the SMI only r#n*ct> variaUon In the toll noiiture condition 
between wilting point end field cepecity. Peck ft el., (19S1) have noted that 
an eapirical relationship between SMI and renotely sensed soil vcitture could 
be developed to sake use of resotely sensed aeasureBents. However, because trve 
SMI only reflects aoitture variations between wilting point and field capacity, 
additional voisture variations in the range of field capacity to saturation 
■easured by reaote sensing could not be reflected in the SMI. 

The use of the i} - tB curves within SSARR would require that the eoisture 
state variable be the actual acisture content ranging froa wilting point to 
saturation rather than the current SMI. One possible approach to revising 
the aodel structure is shown in Fig. 24. The tl • LS curves would define the 
separation of surface runoff, RS, fron the total rainfall infiltration at 

I. 

the soil surface. As long as the average moisture condition was less than 
the field capacity, the rainfall infiltration would be considered to replenish 
soil moisture deficits in the surface zone. Additional infiltrated moisture 
exceeding the soil field capacity would be subiect to deep storage percolation 
(BF) and near-surface lateral flow and inte (RSS) . An additional table 

would define the division of base flow (BF) and interflow (RSS) as a function 
of the soil moisture storage exceeding the soil field capacity. 

In using this approach, the division of rainfall into replenirhed near- 
surface soil moisture, surface runoff. Interflow and baseflow is accomplished 

a 

through two tabular relatii^nships rather thaii the current three* tables. In 
addition, both relationships wt.'.>!d be a function of a single sTate variable. 
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voluMtrlc toil aolkturt cgnttnt rtthtr thtn tht two ytrialUoi SMI and BII. 
The first table (dl - AS curvesl would be defined by the soil texture and the 
current aoisturc conditions at the beginning of the rainfall. The second 
table is Bore eaplrical in nature and would need to be defined through a 
calibration process using rain fall -runoff data froa the particular watershed. 
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THCXJGKTS ON PROBLEMS OF SPATIAL VARIABILITY 

The technique presented in this paper and illustrated in Table 7 estinates 
infiltration capacities at a point. Representative soil hydraulic properties based 
on the SCS soil groups and a reentely sensed estimate of the surface soil moisture 
serve as inputs to the model. As shown in Fig. 2S, taken from a parallel study 
conducted by Nilkening (1981), significant errors can result if the estimate of 
surface soil moisture is in error. Figure 25 is based on a correct surface soil 
moisture of .25 cm/cm. If the hydrologist incorrectly assumes the soil moisture 
to be higher or lower, he would get the indicated error based on experiments \ *.th 
numerical solutions of the Richard's equation. For example, when the moisture 
content is actually .25 and he incorrectly estimates it to be .28 for an analysis 
with the Richard's equation, he would compute the surface runoff to be approximately 
60% higher than actual for a rainfall of 5.1 cm (2.0 inches). Similar patterns 
are encountered for improper estimates of hydraulic conductivity, porosity, etc. 

It is widely recognized that even though a soil type may be consideied to 
be relatively homogeneous, significant variations in soil moisture and hydraulic 
properties over an area can be expected to occur. The sensitivity illustrated 
by Fig. 25 and the recognized t.atura'. spatial variability in soil moisture ana 
soil properties can be quite important in developing a spatially distributed 
hydrologic model. The general approach i<i a hydrologic model is to compute a 
represer;tative infiltration rate and rainfall excess for a discrete areal unit 
rather than for a point. The trends exhibited by Fig. 25 can be tlsed to 
estimate the spatial resolution or cell size to be used in a distributed hydro- 
logic model. 
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Figure 25 . percentase variation in excess raii^all as function of error in 

ESTIfV^TE OF ANTECEDENT fCIS^E 
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Flfure 26 it an idaallxad r«pr«itntation of a watarihad. Two trihutariti 
■orge Into a aain channal draining a «r«'.:arshad having a toil of gtntrally unifora 
toxture. Thare is, however, spatial variability in soil swistura and hydraulic 
properties. Because of variations in thr land use and structure of a aicro- 
drainage network, the hydrologist say want to represent the watershed as an 
an-ay sf cells linked in accordance with Fig. 26. However, if one uses a 
■ean soil aoisture aeasureaent froa soae reaote sensing platfora or a aean 
hydraulic paraaeter based on the soil texture, large errors in the estimate of 
rainfall excess aay result in aany of the cells because of the inherent spatial 
variations of these parameters. As the cell size is increased, however, the 
resulting areal average of the soil aoisture within the cell should approach 
the aean as a representative value. Thus, the desired spatial resolution may have 
to be increased to that shown in Fig. 26b in order to maintain some level of 
precision in the estimate of the rainfall excess. 

This concept of defining errors in terms of spatial resolution has been 
examined, on a preliminary level, using field data representing an example of 
soil aoisture variation within a 23 acre watershed inslruaented by the USDA 
in Chickasha, Oklahoma. The watershed is well-managed pasture with soil variations 
occurring within one soil texture class. Soil aoisture conditions in the top 
five centimeters were represented by the contour map shown in T-ig. 27. 

Using a 17 acre rectangular area within the watershed, a 4Sx48 retangular 
grid matrix was overlaid onto the contour map and encoded cell by cell to 
represent the aoisture variation as an array of 2034 cells. ^ Each cell was 


71 




ALLOM FOR PRECISION IN ESTIIMTION OF 
RAINFALL EXCESS 







SOILS AHO T0P06R;»PHIC wap 


FIGURE 27. WATERSHED R5 


SOIL MOISTURE VARIATION (top 5 cm.) 
MAY 10, 1978 

CHICKASHA, OKLAHOMA (23 ACRES) 
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approiiaatclx 7 ■ x 4 ■ and rapraaantad about 0.04% of tha antira 17 acra araa. 

Tha Man toll Miatura «ai calculatad as .2S cm/cm and tha coafflciant of 
variation was 14% tor this data sat. 

Bacausa the contour Mp onl)r reflects general trends in aoistura variation, 
two additional distributions with coafficlants of variation of 17% and 24% ware 
dailvad frow tha original Mtrlx of calls by allowing each aoistura value to vary 
by introducing a random component through the use of a Markov chain, figure 
28 shows tha symbolic Mps for coefficients of variation of 14* and 17%. 

Assume that the hydrologist wants to represent this watershed as a 
spatially distributed hydrologic model arranged as an array of cells. The rainfall 
excess in each cell will be computed based on the use of the Man soil moisture of 
.25 cm/ cm. However, we do not want the coieputed rainfall excess to deviate from 
the actual theoretical rainfall excess by more than 10% in any cell. If we use 
the three hour storm of Fig. 25. we can accept an error of lO.OS cm/ cm (Range: 
.22-.28cm/ca) in the soil moisture and still esi'lMte a rainfall excess that will 
be within 10% of the correct valve. Thus the question becomes, how largo must our 
spatial resolution element be in order that the average moisture within the 
cell be wihtin 7.03 cm/cm of the Man soil moiiiture. 

By aggregating the cells into new rpitial units such as 2x2, 5x3, 6x6, etc., 
matrices some interesting insights were gained. For the area having a coefficient 
of variation of 14%, if we attempt to model infiltration in each of the 2304 cells 
bated on the assumption that the Man soil moisture of .25 is representative for 
that cell, the computed rainfall excess would be in error by Mre than 10% in 40% 
of the cells. If the spatial resolution were increased to .13 acres C<75% of 
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th* total araa) 33% of tha col la would oacood tho 10% orror lialt. 

If tho roaolutlon mi afaln incraatod to .31 acroa C3% of tho oyorall aroa] 20% 
of tho colls would oxcood tho Halt. Finally, if tho coll rotolution was 
Incroaaod to 1.17 acros (11% of tha total aroa) tha rainfall axctit would bo 
prodlctod within 10% of tho actual values in all of tho colls. In siailar 
ozaainationt for tho coofficlonts of 17% and 24% results were nearly identical. 

It is recofnised that this ozaaination of spatial variability is very 
prelialnary and that sifnificant additional work aust be undertaken. However 
It does illustrate an iaportant consideration that hydrologists are going to 
have to aaka when they consider structuring spatially varied aodels. For 
axaaple, suppose the resolution eleaent of a sensor estiaating surface soil 
aoisture is Halted by technology or econoaics to a .5 kiloaeter square. 

This would give an area of .23 kiloaeters or approxiaately 62 acres. Thus, 
the sensor would give us the aean soil aoisture for a 62 »cre cell. The 
hydrologist could use this aean soil aoisture as the input to the aodel il- 
lustrated in Table 7 to develop the rainfall excess as for a spatially dis- 
tributed aodel having 62 acre cells. However, he would have to look at the 
coefficient of variation and consider the probleas outlined in this chapter 
if it ms necessary to further subdivide the mtersbed into saaller cells to 
siaulate the consequences of spatial variability in other hydrologic processes. 
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conclusions 

Rainfall Infiltralon rataa ara aanaitlva to tha antacadant soil aoiiturc 
conditions and che hydraulic propartias of tha soils. Thaaa sansitivitias 
hava baan invastigacad through tha usa of nuBarical solutions of partial 
diffsrantial aquations that dascriba flow through unsaturatad porous vadia. 
Although thssa nusarical solutions provide good astisatas of tha infiltration 
ratas undar a varlaty of conditions, thair usa in practical hydrologic sK>dels 
dasignad to siaulata watarshed runoff is prohibitad by the coaplaxity of their 
solutions and axtansiva coapjtiv running tlaes. 

A series of functional relationships are presented that allow the 
astiaation of tiae varying infiltration capacities for an array of Initial 
aoisture conditions and soil hydraulic properties. Tha procedure was 
developed using nuaerical axperiaents with the Richard's aquation for 
initial SK>irtura and soil textures t!,at can be expected in watershed hydrology. 
The app' oach is designed to utilise raacww sensing techniques to define the 
iv^itial soil aoisture and Soil Conservation Service soil survey aaps to 
estiaate the soil properties. 

Exaaples are presented that show infiltration ratas predicted by this 
relatively siaple technique coapare very closely with those computed with 
nmarical sensitivity to the Initial soil aoisture canters on that aoisture 
stored in the surface layer, approxlaataly five cantiaaters, of the soil 
coluan. Below that, there is an influence created by variations In the soil 
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Bolitur* prof lit, but tho ijipoct of tho profilo variation li ralntivoly a»all 
by coaparlton vlth othar uncortalntlas that on* would azpoct in tha mtarihad. 

The procedure presented has aeveral key advantages: 1) It is straight 

forward and can be l^>lMsented by the practicing hydrologist using a desk 
top procedure or a very efficient coaputer prograsi; 2] it follows the faailiar 
infiltration- storage concept and, therefore, can be bridged into 
a nuBber of current aodel aethodol ogles; 3) it logically incorporates the 
iapact of initial soil aolsture conditions on the watershed response in terns 
of the physical principles governing the process; 4) the paraaeters are physically 
based and can be defined operationally froa the SCS hydrologic soil group 
classifications presented in soil survey aaps; and S) it is structured to 
utilize near surface soil aolsture that can be developed froa reaote sensing 
techniques. 
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RECQIMENDATIONS 

Th« following roc OMcndat ions might bo considorod for tho furthor develop- 
nt of tho approach takon in this paper: 

1) An attempt to modify tho infiltration cooponents within sovoral 
current models such as IISDAHL and SSARR and then apply these models 
to watershed applications wuuld be worthwhile. Such an effort would 
indicate the amount of work required in replacing the current 
infiltration ccsqionent with the - AS curves. In addition, the 
real benefits could be evaluated in watershed applications. 

2) Finally, an attempt could be made to develop siailar approaches 
for a more complex representation of the soil system Including: 
stratified soil properties within the soil column; vegetative root 
zone and surface landcover effects; the spatial variability of rnois* 
ture conditions and soil properties even with a "homogeneous" 
modeling unit. Attempts could be made to apply the AI - AS curves 
given in this study to more comple - systems. For example, perhaps 

a combination of functions cojld be used to represent a stratified 
soil system. In addition, a more co^>lex theoretical model might 
be used to incorporate the effects of surface vegetation «uch as 
the finite element solution given by Neumann et al., (197S) and 
Feddes et al., (1978). 
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APPENDIX A: 


REPRESENTATIVE AI - AS CURVES 
FOR THE FOLLOWING SOILS: 


Ttxture 
\M»my Sand 
Sandy Loaa 
Loan 

Sandy Clay Loan 
Clay Loan 
Silty Clay 


SCS Soil Croup 
A 
B 
B 
C 
C 
D 



Ali INFILTRATION VOLUME (CM)/15 MINUTE PERIOD 


2.3. 


2 . 2 - 

2 . 1 - 


Lo—y Sind (SCS A) 

■ 0.102 ca/ain • 2.41 in/hr 
porosity (Qj) ■ 0.410 ca/ca 
field capacity f^f^) * 0.12 ca/ca 
wilting pt (*yp) ■ 0*®5 ca/ca 




—I 1 T- 

0.5 1.0 1.5 


-1 1 r— — r- 

2.0 2.5 3.0 3.5 4.0 


aS (cm) 


A I CM /15 MINUTE PERIOD 





AI INFILTRATION VOL (CM)/5 MIN. PERIOD 


Loaay Smnd (SCS A) 



AS (cm) 
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A I CM /15 MINUTE PERIOD 





AI INFILTRATION VOL (CM)/5 MIN. PERIOD 


Loaay Smnd (SCS A) 



AS (cm) 
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M INFILTRATION VOL (cm) /IS MINUTE INCREMENT 



.5 


DEPLETED STORAGE FROH INITI 
(CUNULATIVE INFI 


'--1^. 1H"I 


Soil: SCS Hydrolofic Soil Croup B 

Representative Loaa Soil 

K ■0.01 ca/ain. ■ 0.24 in/hr 
s 

Voluaetric aoisture 

saturation 0.4S ca/ca 
field capacity 0.2S ca/ca 
Milting point O.IS ca/ca 

Initial Unifora Moisture Conditions: 



2.5 3.0 3.5 il.O 


I 





COfIDITIONS ^(Cli) 


A I INFILTRATION VOLUME (CM)/15 MINUTE TIME PERIOD 


0.8 



Of POOR Q'jALirY 




AS (cm) 



I INFILTRATION VOL (CM) /15 MINUTE INCREMENT 



■ 0.004 ca/ain ■ 0.094 ia/kr 
porosity (t^) ■ 0.4S ca/ca 
fiald capacity ca/ca 

viltiif pt. (9^) ■ 0.19 ca/ca 


/VfTKaSfT PblSTlflE 


O 



A I INFILTRATION VOLUME (CM)/15 MINUTE PERIOD 










Al INFILTRATION VOL (CM) /15 MINUTE INCREMENT 



0 
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Infilatration Vol (cm) /5 Minute Increment 
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Ai, Infiltration Vol (cm) / 15 minute Increment 


0.7 


0.7 

0.6 



K - 0.002 ca/ain > 0.047 In/Iir 

s 

porosity (0^) ■0.41 ca/ra 
field capacity ca/ca 

vlltiiif point fO^I ■ 0.26 ca/ca 







A iFILTTVkTION Vx (CM) / 5 fliNJTE I 


0 . 



SILTY CUY (SCS D) 

S Mimit* Increarats 


O O 
^ 3 



ANTECEDENT miSTURE 

0.3^1 

0.40 

0.45 



1.4 


1.6 


1.8 
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OlflCKSlOK *>'01 ST (6>t*C40)«*NCi>C40>»CU»eC40) •T(40),«(40) 
RIAL KSAT 

C0M>*0*/f UNC/NPU tOlNf (6«4'') *0$ TO ( 6| 4C ) tS TOR (4G> 


RIAD SOI 
TMI 
TH§ 
OIRTH 
KSAT 


SOIL ROISTURE RATA 

i K«5| 


UAATEb MOISTURE COUTCKT IROROSITT) C»/CM 
OIPTH Cl AVAILARLI STCRACf INCHES 
SATURATED CONDUCT IVl 1 ; CH/NIN 


READ(S*TCQC) TMliTNSt DEPTH, KSAT 
READ INULTRATION-STORACC RELATIONSHIPS 
NCUR NO. or CURVES PRCVIOED 

NPTS KO. or values OrSCRiaiNK EACH CURVE (TABULAR) 

DELT TIKE InCSEPCNT USED lOR SOLUTION (HOURS) 

RHOIST(I) initial PiOlSTURr CONDITION CP./CP 10* CURVE 1 

D|NM1,J) INriLTRATlON VOLUME PER PERIOD - JTH POINT ON CURVE I 

DSTOdfJ) CHAN6E IN INITIAL STORAGE • JTH POINT ON CURVE I 


READ(S,1CC0) NCUR.N^TS.DELT 

RE AO (Stic 30) (RP0l$T(l),l«1 ,NCUR> 

REAO(S,lLrC) ((DST0(l,j),01Nr(I,J 

L READ RAINIALL CONDITIONS 
C N KUHIER or RAINIALL PERIODS 


J),J*1,NPTS),1»1 ,NCUR) 


RDELT DURATION 01 EACH RAINIALL INCREMENT HOURS 
R(I) RAINIALL RATE lOR ITH INCREMENT INCHES/HOUI 

REA0(S,1CfQ> N,RDELT 
READ(S,10r0) (R(1),I>1,N) 


C convert RAINIALL DATA SO THAT TIME INCREMENTS EQUAL DELT 

II (RDELT.FC.DElT) 60 TO 100 

DI V»kCELT/DELT- 

NDIV>D1V 

RN«N 

RN«01V*PN 

N>RN 

C 

DO 1y 1«1,N 
DO 11 1I«1,ND1V 
RNEU(11)>R(I) 

11 continue 

1C CONTINUE 

C 

DO 12 1«1,N 
R(l)«RNEk(l) 

12 CONTINUE 
C 

ICU DO 14 1>1,N 

R (l)«fi(I )*2.S4*0ELT 
U CONTINUE 

C 

^ CALCULATE MXIMUK AVAILABLEE STORAGE DURING INFILTRATION EVENT 

SMAX»(THS-TH1) •DEPTH*2.S4 
C 

C DETERMIN CURVES TO BE USED FOR INTERPOLATION PROCEDURE 


DO 4S0 I«1,NCUR 

if(thi.gi.rmoist(i>) go to 400 

f lurvNi 

ICURVl«I-1 

R1>I 

IKI.EO.DGO TO 401 

RATIU«(RPOIST( ICURV2)-TH1)/(RM0IST(ICURV2)-RM01S1(ICURV1 )) 

CO TO 45L 

RAT|0>C. 

IF ( I .LT.NCUR) GO TO 4S0 
iCURV2>I 


NCUR) GO TO 4S0 



so. 


4S0 

c 


nV!t!o° 

CONTUUC 


^POOR 




c 

c 

I 

c 


c 

1 

79 


sro 


129. 


calculate AAlltKL CICESS 70R EACH linC INCDIMCNT 
S10A< 

cut 

TC 

t?(S?0Rn?i!6T!srA|) CO TO 502 

CALL SUBAOUTILC TO C/LCULATC INriLTAATlOM CAPACllT f OA 
CAISTING SOU 5 T 0 AAGE 

THl ABTCCCDEnT MOISTUAC CB/CM 

STOR OFPLCTLO STOAAGC FROB INITAU COK&ITIOBS CB 
1 CURV 2 Ft'KCTiOB FOR BOISTURE COBblTlON CLOSEST TO TNI ( > TND 
1 CURV 1 FUNCTION FOR BOISIUME CONDITION CLOSET TO THI ( < THl) 

RATIO INTERPOLATION FACTOR IN USING POTN FUNCTIONS 
DELI INFILTRATION VCLUBE PER T 1 B£ PERIOD CP 

CALL INTERP (TM I . 1 T , 1 CURV2 * 1 CUR Vi , R AT 1 0 tO EL 1 ) 

GO TO 503 
DCLi>KSAT*DELT«60. 

GO TO SOA 

IF (DELI.CT.Pdl-D) DCL 1 >R( 1 T- 1 ) 

ST 0 R( 1 T)>ST 0 R(IT- 1 ) 4 DEL 1 
WRITE ( 7.2 200 ) STORdT > 1 ) ,ST 0 R(IT),DEL 1 
FORBAT ( 319 , 3 ) 

0 ( 1 T- 1 )>(R()T> 1 )-DEL 1 )/?. 3 A 
CUB 0 ( 1 T)>CUB«.( IT -1 ) 4 G( 1 T- 1 ) 

T( 1 T)»T( 1 T- 1 ) 4 DELI 
CONTINUE 
DO 70 C I 5 » 1 ,N 
R( 1 S)«R(IS)/(DILT 42 ,S 4 ) 

STOR(IS)*STOm( is>/ 2 .Sa 
CONT INUE 

RNEWd )>R «1 )>DELT 
DO 701 XS» 2 ,N 

RNEW(IS)>RNCw( IS)«r. ( 1 S)*DELT 
EONTINUE 

WpiTL Mlli (T(I)(R(I),hSEW(I)tSTOR(l>tW(l),CUBG(I),lBl,N> 
FORPAT O 

FORBAT (// 1 k.'T 1 PE(HRS)'. 5 |. **RAINFALL 1 NT ENS I T T < 1 N/ HR ) ' , 

* 51 , 'rainfall VOL(lN)',f S,*INFILTRATIC n VOL(IN) .^At 

• 'rainfall EXCESS! IN)', 5 k. 'CUBULATJVE HUNOF F ( IN) *) 
rOR*'AT( 6 (FlC. 3 , 1 Ck)) 

STOP 
END 


502 

503 
50i 

2000 


500 


700 


701 

C 


icjo 


FTN 2U IPANX 4.72 DBANN 521 COBPON 


i 


^^R V°"oiWi/n-10:5 IMkRPd!,) 


k 

3. 

A. 


SyOROUTlNE iNTERPdHl.lT.lCURyi.lCURV 1,RAT 
CONBON/FUNC/NPTS,DINMu,Zn,DST6(6,A0),S1O 
DIBFNS ION DEL(2) ,IC(2) 


5 

6 
7 
S 
9 

10 

11 

\l 

u 

15 


C 

C 


DO 5000 1>1,2 
IF (IC(I).EO.U) CO TO 5000 

DO 5001 J>2|NPTS 

IF(STCH(1T-1).6T.DST0(IC(I),J)> GO TO 5002 
IF(STOR(lT-1).EQ.DSTO(ICm ,J)) GO TO 5003 
DEL(1)>DINF (1C(I),J-1)«(ST0R(1T-1)-DST0(IC 
« (OlNF(IC(I),J)-DINF(lC(l>,J-1))/(DSTOrjC 




{!> !j) UsTOC 1C (I ) 




1 )) 



1C01 
S0C3 
II 


«0 TO 5CC0* 



RklMfi 16 
fTN 1 WARNINGS 


ir (J.IO.r.R^S) »*LCX)«61*«M| Cfl >,i) 

^S»i!!!j! 

If UCURVI .EQ.O) 60 TO ^00 
ML1*6IL <l) 'RATIO* (Oil Cl >-DEL(1 )> 
60 TO 5110 
OELMDFL Cl) 

fORHATUflO 5)^^*“* 

RETURN 
END 

OUMNT AR6UNCNT 'INI' IS NEVER RCffRENCED 
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